Neutral Plasma Oscillations at Zero Temperature 



(N 
O 

o 

u 

Oh- 
< 

(N 



c*3 

a 

11; 

1/3 

o 

C/5 



> 

00 

o 
o 
o 

1/3 

o 

<Z3 1 

^ : 

OH- 



X 



S. D. Bergeson and R. L. Spencer 
Department of Physics and Astronomy, Brigham Young University, Provo, UT 84602 

(Dated: February 2, 2008) 

We use cold plasma theory to calculate the response of an ultracold neutral plasma to an applied 
rf field. The free oscillation of the system has a continuous spectrum and an associated damped 
quasimode. We show that this quasimode dominates the driven response. We use this model to 
simulate plasma oscillations in an expanding ultracold neutral plasma, providing insights into the 
assumptions used to interpret experimental data [Phys. Rev. Lett. 85, 318 (2000)]. 

PACS numbers: 52.55.Dy, 32.80.Pj, 52.27.Gr, 52.35.Fp 



Recent experimental [i] p, |], ||, |(| and theoretical 
work H g, [lO], [n], [12, fUjl has studied the formation 
and evolution of ultracold plasmas. In the laboratory, 
cold plasmas are created either by directly photo-ionizing 
laser-cooled atoms, or by exciting the atoms to high-lying 
Rydberg states that spontaneously ionize. A fraction of 
the electrons escape the plasma and the resulting electric 
field drives the ion expansion. Models of the expansion 
suggest that the density profile in the plasma is approx- 
imately Gaussian, and that it expands in a self-similar 
manner. It can be expressed as 



n(r,t) 



exp[ 



2 /2a 2 



(1) 



where a = \J u 2 + v 2 t 2 is the time-dependent width of 
the distribution, v is the asymptotic expansion velocity, 
and t is time. 

The experimental justification of this density profile for 
the case of an expanding plasma is based on the plasma's 
response to a spatially uniform applied rf field. In 
those experiments, Xe atoms initially cooled to ~ 
were ionized by a dye laser. The initial electron energy 
(E e /ks) ranged from a few to 1000 K, and the initial 
electron density ranged from 0.2 to 2.5 x 10 9 cm~ 3 . The 
plasmas were nearly charge-neutral, and at electron en- 
ergies above 70 K, the kinetic energy of the electrons al- 
lowed significant loss. The resulting net positive charge 
in the cloud drove the plasma expansion. 

As the plasma expanded and its density decreased, the 
applied rf field pumped energy into the plasma. The 
heating was assumed to be largest where the applied rf 
frequency matched the local plasma frequency. Because 
of collisions in the plasma, the local heating presumably 
raised the overall plasma temperature, and a small num- 
ber of the more weakly-bound electrons were ejected. 
The experiment measured the rate at which electrons 
were ejected from the expanding plasma as a function 
of time for a fixed applied rf frequency. This signal was 
presumed to be proportional to the rate at which the rf 
field heats the plasma, and was a measure of the weighted 
time-dependent density profile. 

The peak of this signal was interpreted to correspond 



to the "average" density of the plasma, 



(t) = 2- 3/2 n(0,i) = 



n crl 



3/2 



(2) 



with v a constant. This density was experimentally de- 
termined by setting the applied frequency oj equal to 
the average plasma frequency, u> p , and using the rela- 
tion uj p = yj 'q 2 fi(t) /m e €o, where q is the electron charge, 
m e is the electron mass, and eo is the permittivity of free 
space. The derived density, n, is time dependent. For 
a different applied rf frequency to, the signal peaks at a 
different time. A least-squares fit of fi{t) to Eq. (J|) gives 
the initial density no and the expansion velocity v. 

In this paper, we examine the validity of the assump- 
tions used in this interpretation by presenting the re- 
sponse of the system predicted by cold plasma theory. 
This theory is good for plasmas in which (\r>/a) 2 <C 1, 
where Ad is the Debye length. In the experiments we are 
simulating, (\r>/<j) 2 ~ 10~ 3 — 10 _1 , making cold plasma 
theory a good way to model the collective response of the 
system. 

We consider spherically symmetric charge-neutral 
Gaussian distributions of cold ions and electrons to which 
an rf-electric field in the z-direction is applied. The ions 
are taken as fixed on the time-scale of the rf oscillation. 
The fluid equations describing the electrons are 



On 
~dt 



+ V • {nv) = 



V <j> = -q{n - rii) /e 



dv 



m | — — h v ■ Vu = qE — jmv. 
ot ' 



(3) 
(4) 
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In these equations n = n{r, t) is the electron density, 
rii = rii(r) is the ion density (assumed to be fixed in 
time), v — v(r,t) is the electron velocity, E = E(f,t) 
is the total electric field (the sum of the external field 
and that generated by the plasma response), and 7 is 
a phenomcnological damping rate. Some effects of fi- 
nite temperature, such as electron-electron collisions and 
electron-ion collisions, are approximated by the the pa- 
rameter 7. 
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FIG. 1: Quasimode frequency, ui m (o), and damping rate, 
7m (o), for density profiles given by Eq. (^). For p = 2, 
the distribution is Gaussian. As p increases, the distribution 
approaches a "top-hat," the quasimode frequency increases, 
and the quasimode damping decreases (see text). 



In the experiment by Kulin et al. p{, the applied E 
was uniform in space and oscillating in the ^-direction. 
Assuming that this applied field is small, the fluid equa- 
tions can be linearized by assuming that the density 
and the potential are of the form n(r, 9, t) = (r) + 
Sn(r) cos 9e~ luJt , where n^(r) is the equilibrium elec- 
tron density, and 4>(r, 9, t) = 5<j>(r) cos 9e~ tuJt (spherical 
harmonics with I = 1). After some algebra to linearize 
the fluid Equations the potential produced by the 

electrons is given by the expression 



C64> 



lu 2 — ujp(r) + iuj^i 



2 — ujp(r) + iujj 



where the operator C is defined as 



dr 2 



2 d 
r dr 
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and where w p (r) = y q 2 n^ (r) /m e eo is the plasma fre- 
quency. The boundary conditions are that 6<j>(0) = and 
that at infinity the field is that of a dipole, 6<j)(r) oc 1/r 2 . 

Because this is a driven system, we first look for nor- 
mal modes. Notice, however, that with Eq = and 7 = 
(undamped free response) this mode equation has a con- 
tinuous spectrum, similar to the diocotron mode equation 
in non-neutral plasmas JTJ, ^ Q [fij [D| And, as 
with the diocotron mode, this equation also has a damped 
quasimode. Following Ref. [l5| , Eq. (||) can be solved 
along a contour in the complex r-plane to uncover this 
damped quasimode. Figure shows the frequency and 
damping rate of this quasimode for a sequence of density 
profiles of the form 



FIG. 2: Power absorbed by the (non-expanding) plasma as a 
function applied frequency. The solid line shows the absorbed 
power for a Gaussian density profile, p = 2 in Eq. (^). Also 
shown is the absorbed power for super-Gaussian profiles, with 
p = 4, 6. The mode damping decreases with increasing p. 
These plots are generated with 7/ojp ' = 0.01. For larger 
values of 7, these curves are somewhat broader and shifted to 
slightly higher ui. 



The real and imaginary parts of the quasimode fre- 
quency are well approximated by the simple formulas 
Lu m /uj p (0) = l/\/3-0.668/p-0.102/p 2 and 7 m /w p (0) = 
1.077/p- 0.959/p 2 . 

Note that the Gaussian (p = 2) quasimode is heavily 
damped and that as p becomes large (making the density 
profile approach a step-profile) the damping rate goes to 
zero and the frequency approaches u m = u) p (0)/V3. 

The quasimode dominates the free response of a sys- 
tem with a continuous spectrum [^0| |2l[] . It also strongly 
influences the driven response, as can be seen in Fig. |^. 
This figure shows the rate at which a plasma absorbs 
power from an applied rf field as a function of the fre- 
quency of that field. To calculate this power, we note that 
the total electric field (the applied field plus the plasma 
response) is given by 



E(r, t) = (zE cos 9 - V50) e" 



(9) 



By using the linearized Eq. (Q) to compute the current 
density J = qnv in the plasma, the power density can be 
written as 



(P) = \Re 



E-.l I = ^ 



E 



uj 2 + 7 2 



(10) 



n(r) = n e- rP ^ P 



(8) 



The influence of the quasimode on the absorbed power 
can be seen by integrating this power density over all 
space for various values of the driving frequency u>. Al- 
though strongly damped, the quasimode still influences 
how the plasma responds to the rf field. Figure |^ shows 
the calculated power absorbed by the plasma as a func- 
tion of the applied rf frequency for a non-expanding 
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plasma. In each case the peak of the absorbed power 
occurs near the frequency of the quasimode. Notice that 
these resonances occur at quite small values of the den- 
sity, rather than in the body of the density distribu- 
tion as assumed by Kulin, et al. The experiment used 
Lu/uj p (0) = 0.6 

We now use these cold driven solutions of Eq. (||) to 
simulate the plasma expansion experiment. In order to 
do this, we need an expression for the electron density, 
n(r,t) and the phenomcnological damping rate 7. For 
both of these we use the model of Robicheaux and Han- 
son (ll[] . They simulated an expanding ultracold neutral 
plasma to study three-body recombination and electron 
heating. An important result of their paper is that the 
electron temperature scales with density as T e cx n~ 1//3 , 
which makes the Coulomb logarithm (In A) constant. We 
use this result in our simulation of the collective response 
of the plasma to an applied field. It is likely that the phe- 
nomenological damping rate, 7, is proportional to the 
electron collision rate, v. The scaled collision rate in the 
plasma is vjio v cx -y/n/T 3 In A, which is constant in this 
temperature scaling. For typical experimental param- 
eters, this gives j/oj p (0) w 1/100. Another important 
result from the work of Robicheaux and Hanson is that 
most of the ions in the plasma experience an electric field 
linearly proportional to r. Given the known initial condi- 
tions in the plasma, this leads directly to the result that 
the plasma density is Gaussian, and that it expands in a 
self-similar manner described by Eq. (|l|) |2^| . 

We simulate the plasma expansion experiment as fol- 
lows. The plasma expansion velocity (v), the rf frequency 
(u>) and rf amplitude (Eo) are held constant. At a partic- 
ular time t, we insert the density from Eq. (|l|) into Eq. 
(^) to find the potential 5<p and the total electric field 
[Eq. (0)]. We integrate the power density [Eq. fli~0[)] over 
the spatial coordinates to get the plasma heating rate 
at that time. We increment time and repeat the heating 
rate calculation to generate one of the time-sweeps shown 
in Fig. ||(a). Repeating these calculations for a range of 
applied field frequencies generates all of the time sweeps 
shown in the figure. 

We convert the time sweeps in Fig. ||(a) into aver- 
age density determinations. At the peak of each time 
sweep we set the applied frequency lo equal to the av- 
erage plasma frequency uj v — y^q 2 h/m e eo to calculate 
the average density n, as was done in the experiment 
With this conversion from frequency to density, the 
simulated average density n(t) is plotted in Fig. ||(b). 
The data is fit to Eq. (g) using a least-squares method 
to extract the initial central density, no, and expansion 
velocity v. 

The shapes of the calculated plasma heating rates 
shown in Fig. ^(a) match the experimental ones from 
Reference B], and the fit of n(t) in Fig. ||(b) to Eq. (||) is 
excellent p3j . The expansion velocity extracted from the 
fit exactly matches the known velocity that we put into 
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FIG. 3: (a) Calculated plasma heating rate. The applied fre- 
quencies ui/2n in MHz are shown to the left of each trace. 
The traces are offset vertically for clarity. The plasma con- 
ditions are no = 10 15 m -3 , v = 100 m/s, 00 = 220 um, and 
7/a>p(0) = 0.01. (b) Average plasma density as a function 
of time. The circles (•) are taken from the data in (a). The 
applied frequency is converted to a density, n(t), with t given 
by the peak of the signal (see text). The solid line is a fit of 
these points to Eq. (2h. The fitted density is n = 1.1 x 10 14 
m -3 or no = 3.1 x 10 m -3 , about a factor of 3 too low. The 
fitted expansion velocity is v = 100 m/s. 

the simulation, independent of the value of the scaled 
damping rate "f/u) p (0). If the Robicheaux and Hanson 
temperature scaling is correct, the plasma response to 
the rf field accurately reproduces the true asymptotic ex- 
pansion velocity. However, different temperature scalings 
make the fit worse. If instead of using T cx n 1 / 3 in deter- 
mining 7, we use T cx n 2 / 3 , the fitted velocity increases 
by roughly a factor of two, depends on j/ui p (0), the ex- 
tracted n(t) does not follow Eq. (0), and the simulation 
does not fit the experiment. 

The initial central density no extracted from the fit of 
n(t) to Eq. (0) is quite sensitive to the scaled damping- 
factor, j/u} p (0). It scales approximately as e 27 /" p< - - ) . It 
also depends on the particular temperature scaling used 
in the simulation, and can vary by orders of magnitude. 
For the Robicheaux and Hanson temperature scaling and 
with 7/w p (0) = 0.01, the fitted central density underes- 
timates the known central density in our simulation by a 
factor of three. 

Cold plasma theory only addresses the collective be- 
havior of the system. Another effect that might be 
involved in the resonant behavior seen in the experi- 
ment is that the electron motion in the plasma resonates 
with the driving field. Robicheaux and Hanson estimate 
that in the bulk of the plasma the potential satisifes 
q^/ksT = r 2 /2a 2 . This means that most of the particles 
are in an isotropic harmonic oscillator well with natural 
frequency ojo = u) p (0)\d(0)/ct, where Xd(0) is the Debye 
length at the center of the cloud. If the peaks in the 
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experiment were due to particle resonance at u> = u>o in 
the expanding plasma, the experiment does not indicate 
a uniform expansion p4| . 

However, in experiments where Xd(0) <C <j, such as in 
the present simulations, we would expect collective effects 
like the quasimode, to dominate over individual particle 
effects like a resonance at w = ujq. But even though 
collective effects should dominate these plasmas, particle 
resonance can play a role, as it does in the case of electron 
Bernstein wave in magnetized plasmas p5[ . Calculating 
this effect requires a kinetic theory treatment that is be- 
yond the scope of this paper. But we note that even in 
kinetic theory the particle resonance will not be sharp 
because the well is not perfectly parabolic, causing the 
orbit frequency to vary with the energy and angular mo- 
mentum. Averaging this resonance over the distribution 
function blurs its effect and makes the resonances at u>o, 
2u>o, etc. less pronounced. 

It would be very helpful to have direct measurements 
of the ion density distribution, such as with optical de- 
tection. These measurements could provide benchmark 
data against which future simulations can be compared. 
Also, as mentioned previously, full kinetic simulations are 
probably required to address the particle orbit problem. 
Such measurements and simulations are currently under- 
way in our laboratory. 

In conclusion, we have used a simulation based on cold 
plasma theory to calculate the response of an ultracold 
plasma to an applied rf field. When T oc n 1 / 3 , the sim- 
ulation accurately reproduces the experimental data of 
Kulin et al. providing an important consistency check 
of the model of Robicheaux and Hanson |llj . 

The model confirms that the plasma response accu- 
rately reproduces the asymptotic velocity of an expand- 
ing plasma, and determines the initial plasma density to 
within a factor of two or three. The plasma response is 
dominated by the quasi-mode. Even though the mode is 
strongly damped, it shifts the plasma resonance to lower 
densities than otherwise expected from a local-density 
approximation. 
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